Genomic analysis of male puberty timing highlights shared genetic basis with hair colour and lifespan

The timing of puberty is highly variable and is associated with long-term health outcomes. To date, understanding of the genetic control of puberty timing is based largely on studies in women. Here, we report a multi-trait genome-wide association study for male puberty timing with an effective sample size of 205,354 men. We find moderately strong genomic correlation in puberty timing between sexes (rg = 0.68) and identify 76 independent signals for male puberty timing. Implicated mechanisms include an unexpected link between puberty timing and natural hair colour, possibly reflecting common effects of pituitary hormones on puberty and pigmentation. Earlier male puberty timing is genetically correlated with several adverse health outcomes and Mendelian randomization analyses show a genetic association between male puberty timing and shorter lifespan. These findings highlight the relationships between puberty timing and health outcomes, and demonstrate the value of genetic studies of puberty timing in both sexes

Genome-wide association study with 1000 genomes imputation identifies signals for nine sex hormone-related phenotypes.

PublishedJournal ArticleResearch Support, Non-U.S. Gov'tThis is the final version of the article. Available from Nature Publishing Group via the DOI in this record.Genetic factors contribute strongly to sex hormone levels, yet knowledge of the regulatory mechanisms remains incomplete. Genome-wide association studies (GWAS) have identified only a small number of loci associated with sex hormone levels, with several reproductive hormones yet to be assessed. The aim of the study was to identify novel genetic variants contributing to the regulation of sex hormones. We performed GWAS using genotypes imputed from the 1000 Genomes reference panel. The study used genotype and phenotype data from a UK twin register. We included 2913 individuals (up to 294 males) from the Twins UK study, excluding individuals receiving hormone treatment. Phenotypes were standardised for age, sex, BMI, stage of menstrual cycle and menopausal status. We tested 7,879,351 autosomal SNPs for association with levels of dehydroepiandrosterone sulphate (DHEAS), oestradiol, free androgen index (FAI), follicle-stimulating hormone (FSH), luteinizing hormone (LH), prolactin, progesterone, sex hormone-binding globulin and testosterone. Eight independent genetic variants reached genome-wide significance (P<5 × 10(-8)), with minor allele frequencies of 1.3-23.9%. Novel signals included variants for progesterone (P=7.68 × 10(-12)), oestradiol (P=1.63 × 10(-8)) and FAI (P=1.50 × 10(-8)). A genetic variant near the FSHB gene was identified which influenced both FSH (P=1.74 × 10(-8)) and LH (P=3.94 × 10(-9)) levels. A separate locus on chromosome 7 was associated with both DHEAS (P=1.82 × 10(-14)) and progesterone (P=6.09 × 10(-14)). This study highlights loci that are relevant to reproductive function and suggests overlap in the genetic basis of hormone regulation.We thank Roche Diagnostics Australia Pty Limited, Castle Hill, Australia, who provided support for the analysis of the hormones. We thank the volunteer twins for their participation in the study. Twins UK received funding support from NIHR Biomedical Research Centre (grant to Guys’ and St Thomas’ Hospitals and King’s College London); the Chronic Disease Research Foundation; Canadian Institutes of Health Research, the Canadian Foundation for Innovation, the Fonds de la Recherche en Santé Québec, The Lady Davis Institute, the Jewish General Hospital and Ministère du Développement économique, de l'Innovation et de l'Exportation du Quebec. The Australian National Health and Medical Research Council (NHMRC project grants 1010494, 1048216), and Sir Charles Gairdner Hospital Research (grant PP2009/028). This work was supported by funding from the Wellcome Trust (092447/Z/10/Z) and Medical Research Council (MC_U106179472)

Genome-wide association study meta-analysis for quantitative ultrasound parameters of bone identifies five novel loci for broadband ultrasound attenuation.

Osteoporosis is a common and debilitating bone disease that is characterised by low bone mineral density, typically assessed using dual-energy X-ray absorptiometry. Quantitative ultrasound (QUS), commonly utilising the two parameters velocity of sound (VOS) and broadband ultrasound attenuation (BUA), is an alternative technology used to assess bone properties at peripheral skeletal sites. The genetic influence on the bone qualities assessed by QUS remains an under-studied area. We performed a comprehensive GWAS including low-frequency variants (MAF ≥0.005) for BUA and VOS using a discovery population of individuals with whole-genome sequence (WGS) data from the UK10K project (n=1,268). These results were then meta-analysed with those from two deeply imputed GWAS replication cohorts (n=1,610 and 13,749). In the gender-combined analysis, we identified eight loci associated with BUA and five with VOS at the genome-wide significance level, including three novel loci for BUA at 8p23.1 (PPP1R3B), 11q23.1 (LOC387810) and 22q11.21 (SEPT5) (P = 2.4 × 10-8-1.6 × 10-9). Gene-based association testing in the gender-combined dataset revealed eight loci associated with BUA and seven with VOS at the genome-wide significance level, with one novel locus for BUA at FAM167A (8p23.1) (P = 1.4 × 10-6). An additional novel locus for BUA was seen in the male-specific analysis at DEFB103B (8p23.1) (P = 1.8 × 10-6). Fracture analysis revealed significant associations between variation at the WNT16 and RSPO3 loci and fracture risk (P = 0.004 and 4.0 × 10-4 respectively). In conclusion, by performing a large GWAS meta-analysis for QUS parameters of bone using a combination of WGS and deeply imputed genotype data, we have identified five novel genetic loci associated with BUA

Genomic analysis of male puberty timing highlights shared genetic basis with hair colour and lifespan

The timing of puberty is highly variable and is associated with long-term health outcomes. To date, understanding of the genetic control of puberty timing is based largely on studies in women. Here, we report a multi-trait genome-wide association study for male puberty timing with an effective sample size of 205,354 men. We find moderately strong genomic correlation in puberty timing between sexes (rg = 0.68) and identify 76 independent signals for male puberty timing. Implicated mechanisms include an unexpected link between puberty timing and natural hair colour, possibly reflecting common effects of pituitary hormones on puberty and pigmentation. Earlier male puberty timing is genetically correlated with several adverse health outcomes and Mendelian randomization analyses show a genetic association between male puberty timing and shorter lifespan. These findings highlight the relationships between puberty timing and health outcomes, and demonstrate the value of genetic studies of puberty timing in both sexes

Genetic Predisposition to an Impaired Metabolism of the Branched-Chain Amino Acids and Risk of Type 2 Diabetes: A Mendelian Randomisation Analysis

$\textbf{BACKGROUND}$: Higher circulating levels of the branched-chain amino acids (BCAAs; i.e., isoleucine, leucine, and valine) are strongly associated with higher type 2 diabetes risk, but it is not known whether this association is causal. We undertook large-scale human genetic analyses to address this question. $\textbf{METHODS AND FINDINGS}$: Genome-wide studies of BCAA levels in 16,596 individuals revealed five genomic regions associated at genome-wide levels of significance (p < 5 × 10-8). The strongest signal was 21 kb upstream of the PPM1K gene (beta in standard deviations [SDs] of leucine per allele = 0.08, p = 3.9 × 10-25), encoding an activator of the mitochondrial branched-chain alpha-ketoacid dehydrogenase (BCKD) responsible for the rate-limiting step in BCAA catabolism. In another analysis, in up to 47,877 cases of type 2 diabetes and 267,694 controls, a genetically predicted difference of 1 SD in amino acid level was associated with an odds ratio for type 2 diabetes of 1.44 (95% CI 1.26-1.65, p = 9.5 × 10-8) for isoleucine, 1.85 (95% CI 1.41-2.42, p = 7.3 × 10-6) for leucine, and 1.54 (95% CI 1.28-1.84, p = 4.2 × 10-6) for valine. Estimates were highly consistent with those from prospective observational studies of the association between BCAA levels and incident type 2 diabetes in a meta-analysis of 1,992 cases and 4,319 non-cases. Metabolome-wide association analyses of BCAA-raising alleles revealed high specificity to the BCAA pathway and an accumulation of metabolites upstream of branched-chain alpha-ketoacid oxidation, consistent with reduced BCKD activity. Limitations of this study are that, while the association of genetic variants appeared highly specific, the possibility of pleiotropic associations cannot be entirely excluded. Similar to other complex phenotypes, genetic scores used in the study captured a limited proportion of the heritability in BCAA levels. Therefore, it is possible that only some of the mechanisms that increase BCAA levels or affect BCAA metabolism are implicated in type 2 diabetes. $\textbf{CONCLUSIONS}$: Evidence from this large-scale human genetic and metabolomic study is consistent with a causal role of BCAA metabolism in the aetiology of type 2 diabetes.MRC Epidemiology Unit, Fenland study, EPIC-InterAct study, EPIC-Norfolk case-cohort study funding: this study was funded by the United Kingdom’s Medical Research Council through grants MC_UU_12015/1, MC_UU_12015/5, MC_PC_13046, MC_PC_13048 and MR/L00002/1. We acknowledge support from the National Institute for Health Research Biomedical Research Centre. The research leading to these results has received support from the Innovative Medicines Initiative Joint Undertaking under EMIF grant agreement number 115372, resources of which are composed of financial contribution from the European Union's Seventh Framework Programme (FP7/2007-2013) and EFPIA companies’ in kind contribution. EPIC-InterAct Study funding: funding for the InterAct project was provided by the EU FP6 programme (grant number LSHM_CT_2006_037197). MRC Human Nutrition Research funding: This research was supported by the Medical Research Council (MC_UP_A090_1006) and Cambridge Lipidomics Biomarker Research Initiative (G0800783). The SABRE study was funded at baseline by the UK Medical Research Council, Diabetes UK and the British Heart Foundation and at follow-up by a programme grant from the Wellcome Trust (WT082464) and British Heart Foundation (SP/07/001/23603); Diabetes UK funded the metabolomics analyses (13/0004774). RJOS, EN, JRZ and AK received funding from the Swedish Research Council, Stockholm County Council, Novo Nordisk Foundation and Diabetes Wellness. DBS is supported by the Wellcome Trust grant number 107064. MIM is a Wellcome Trust Senior Investigator and is supported by the following grants from the Wellcome Trust: 090532 and 098381. IB is supported by the Wellcome Trust grant WT098051

Genetic predisposition to an impaired metabolism of the branched-chain amino acids and risk of type 2 diabetes: a mendelian randomisation analysis

BACKGROUND: Higher circulating levels of the branched-chain amino acids (BCAAs; i.e., isoleucine, leucine, and valine) are strongly associated with higher type 2 diabetes risk, but it is not known whether this association is causal. We undertook large-scale human genetic analyses to address this question. METHODS AND FINDINGS: Genome-wide studies of BCAA levels in 16,596 individuals revealed five genomic regions associated at genome-wide levels of significance (p < 5 × 10-8). The strongest signal was 21 kb upstream of the PPM1K gene (beta in standard deviations [SDs] of leucine per allele = 0.08, p = 3.9 × 10-25), encoding an activator of the mitochondrial branched-chain alpha-ketoacid dehydrogenase (BCKD) responsible for the rate-limiting step in BCAA catabolism. In another analysis, in up to 47,877 cases of type 2 diabetes and 267,694 controls, a genetically predicted difference of 1 SD in amino acid level was associated with an odds ratio for type 2 diabetes of 1.44 (95% CI 1.26-1.65, p = 9.5 × 10-8) for isoleucine, 1.85 (95% CI 1.41-2.42, p = 7.3 × 10-6) for leucine, and 1.54 (95% CI 1.28-1.84, p = 4.2 × 10-6) for valine. Estimates were highly consistent with those from prospective observational studies of the association between BCAA levels and incident type 2 diabetes in a meta-analysis of 1,992 cases and 4,319 non-cases. Metabolome-wide association analyses of BCAA-raising alleles revealed high specificity to the BCAA pathway and an accumulation of metabolites upstream of branched-chain alpha-ketoacid oxidation, consistent with reduced BCKD activity. Limitations of this study are that, while the association of genetic variants appeared highly specific, the possibility of pleiotropic associations cannot be entirely excluded. Similar to other complex phenotypes, genetic scores used in the study captured a limited proportion of the heritability in BCAA levels. Therefore, it is possible that only some of the mechanisms that increase BCAA levels or affect BCAA metabolism are implicated in type 2 diabetes. CONCLUSIONS: Evidence from this large-scale human genetic and metabolomic study is consistent with a causal role of BCAA metabolism in the aetiology of type 2 diabetes

Evolutionary Consequences of Altered Atmospheric Oxygen in Drosophila melanogaster

Twelve replicate populations of Drosophila melanogaster, all derived from a common ancestor, were independently evolved for 34+ generations in one of three treatment environments of varying PO2: hypoxia (5.0–10.1 kPa), normoxia (21.3 kPa), and hyperoxia (40.5 kPa). Several traits related to whole animal performance and metabolism were assayed at various stages via “common garden” and reciprocal transplant assays to directly compare evolved and acclimatory differences among treatments. Results clearly demonstrate the evolution of a greater tolerance to acute hypoxia in the hypoxia-evolved populations, consistent with adaptation to this environment. Greater hypoxia tolerance was associated with an increase in citrate synthase activity in fly homogenate when compared to normoxic (control) populations, suggesting an increase in mitochondrial volume density in these populations. In contrast, no direct evidence of increased performance of the hyperoxia-evolved populations was detected, although a significant decrease in the tolerance of these populations to acute hypoxia suggests a cost to adaptation to hyperoxia. Hyperoxia-evolved populations had lower productivity overall (i.e., across treatment environments) and there was no evidence that hypoxia or hyperoxia-evolved populations had greatest productivity or longevity in their respective treatment environments, suggesting that these assays failed to capture the components of fitness relevant to adaptation

Physical and neurobehavioral determinants of reproductive onset and success.

The ages of puberty, first sexual intercourse and first birth signify the onset of reproductive ability, behavior and success, respectively. In a genome-wide association study of 125,667 UK Biobank participants, we identify 38 loci associated (P < 5 × 10(-8)) with age at first sexual intercourse. These findings were taken forward in 241,910 men and women from Iceland and 20,187 women from the Women's Genome Health Study. Several of the identified loci also exhibit associations (P < 5 × 10(-8)) with other reproductive and behavioral traits, including age at first birth (variants in or near ESR1 and RBM6-SEMA3F), number of children (CADM2 and ESR1), irritable temperament (MSRA) and risk-taking propensity (CADM2). Mendelian randomization analyses infer causal influences of earlier puberty timing on earlier first sexual intercourse, earlier first birth and lower educational attainment. In turn, likely causal consequences of earlier first sexual intercourse include reproductive, educational, psychiatric and cardiometabolic outcomes.This research has been conducted using the UK Biobank Resource. This work was supported by the Medical Research Council (Unit Programme numbers MC_UU_12015/1 and MC_UU_12015/2)

Prevalence of deleterious variants in MC3R in patients with constitutional delay of growth and puberty.

CONTEXT: The melanocortin 3 receptor (MC3R) has recently emerged as a critical regulator of pubertal timing, linear growth and the acquisition of lean mass in humans and mice. In population-based studies, heterozygous carriers of deleterious variants in MC3R report a later onset of puberty than non-carriers. However, the frequency of such variants in patients who present with clinical disorders of pubertal development is currently unknown. OBJECTIVE: To determine whether deleterious MC3R variants are more frequently found in patients clinically presenting with constitutional delay of growth and puberty (CDGP) or normosmic idiopathic hypogonadotropic hypogonadism (nIHH). DESIGN, SETTING AND PARTICIPANTS: We examined the sequence of MC3R in 362 adolescents with a clinical diagnosis of CDGP and 657 patients with nIHH, experimentally characterised the signalling properties of all non-synonymous variants found and compared their frequency to that in 5774 controls from a population-based cohort. Additionally, we established the relative frequency of predicted deleterious variants in individuals with self-reported delayed vs normally timed menarche/voice breaking in the UK Biobank cohort. RESULTS: MC3R loss-of-function variants were infrequent but overrepresented in patients with CDGP (8/362 (2.2%), OR = 4.17, p = 0.001). There was no strong evidence of overrepresentation in patients with nIHH (4/657 (0.6%), OR = 1.15, p = 0.779). In 246,328 women from UK Biobank, predicted deleterious variants were more frequently found in those self-reporting delayed (≥16 years) vs normal age at menarche (OR = 1.66, p = 3.90E-07). CONCLUSIONS: We have found evidence that functionally damaging variants in MC3R are overrepresented in individuals with CDGP but are not a common cause of this phenotype

Penetrance of pathogenic genetic variants associated with premature ovarian insufficiency

This is the author accepted manuscript. The final version is available from Nature Research via the DOI in this recordData availability: We used publicly available individual-level genotype and phenotype data from the UK Biobank (https://biobank.ndph.ox.ac.uk/showcase/). Access to these data needs to be requested from the UK Biobank (https://www.ukbiobank.ac.uk/enable-your-research/apply-for-access). This research was conducted using the UK Biobank Resource under application 9905 (University of Cambridge) and 9072 and 871 (University of Exeter). POI genes were identified from the GEL Panel App (https://panelapp.genomicsengland.co.uk/panels/155/). Genes were annotated with pathogenicity constraint metrics from gnomAD v.2.1.1.Premature ovarian insufficiency (POI) affects 1% of women and is a leading cause of infertility. It is often considered to be a monogenic disorder, with pathogenic variants in ~100 genes described in the literature. We sought to systematically evaluate the penetrance of variants in these genes using exome sequence data in 104,733 women from the UK Biobank, 2,231 (1.14%) of whom reported at natural menopause under the age of 40 years. We found limited evidence to support any previously reported autosomal dominant effect. For nearly all heterozygous effects on previously reported POI genes, we ruled out even modest penetrance, with 99.9% (13,699 out of 13,708) of all protein-truncating variants found in reproductively healthy women. We found evidence of haploinsufficiency effects in several genes, including TWNK (1.54 years earlier menopause, P = 1.59 × 10-6) and SOHLH2 (3.48 years earlier menopause, P = 1.03 × 10-4). Collectively, our results suggest that, for the vast majority of women, POI is not caused by autosomal dominant variants either in genes previously reported or currently evaluated in clinical diagnostic panels. Our findings, plus previous studies, suggest that most POI cases are likely oligogenic or polygenic in nature, which has important implications for future clinical genetic studies, and genetic counseling for families affected by POI.Medical Research Council (MRC)National Institute for Health and Care Research (NIHR)QUEX InstituteUniversity of CambridgeCancer Research UKInnovative Medicines Initiative 2 Joint UndertakingNational Health and Medical Research CouncilEuropean Research Council (ERC)Novo Nordisk FoundationIndependent Research Foundation DenmarkDanish National Research Foundation Centr